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I. INTRODUCTION 
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High temperature QCD is the subject of a long-running experimental program of heavy-ion collisions in several 
CN| ■ experimental facilities, including the Relativistic Heavy Ion Collider (RHIC) at BNL, the Large Hadron Collider 
(LHC) at CERN, and planned programs in GSI and Dubna. They aim to collide together heavy-ions at high energies 
and produce fireballs of strongly interacting matter, whose study may then yield information about the phases and 
properties of such matter. Lattice computations are the primary source of theoretical information on such matter. 
The aim of this study is to examine the response of this matter to quasi-static perturbations. 

It is known that the generic response of matter to such perturbations is to screen them. The knowledge of screening 
lengths, £, is one of the most basic pieces of microscopic information we can have about the system. Of particular 
interest is the longest screening length, £o, or its inverse, the smallest screening mass, [j,q = l/£n Q- When = 0, the 
'— effective long-distance theory of matter (at, or close to, equilibrium) requires us to take into account these unscreened 
perturbations. Such is the case in a QED plasma, where magnetic fields are not screened. Not only are equilibrium 
properties of a plasma strongly influenced by this, but also the off-equilibrium long-distance theory changes to magneto 
ly-j , hydrodynamics rather than Navier-Stokes hydrodynamics, as it is for other fluids. 

1 For the QCD plasma with non-vanishing quark masses, one knows that all fields arc screened and the long-distance 
theory will be hydrodynamics coupled to the diffusion of conserved charges. Nevertheless, the study of screening 
■ masses is of practical importance. If the smallest dimension of the fireball produced in heavy-ion collisions is £, 
| one expects thermodynamic properties to manifest themselves only when ifig 1. Furthermore, since both these 
£3 . quantities are functions of the temperature T, it is possible for the fireball to drop into, and out of, equilibrium at 
different temperatures. In this way, quasi-static properties such as the screening masses may put bounds on truly 
i dynamical quantities such as the thermalization and freezeout times in heavy-ion collisions. 
J> ■ The study of screening properties in a plasma has a long history. In the glue sector, the Debye screening length 
has been the object of many studies and now seems to be quantitatively understood, both in non-perturbative lattice 
studies and in weak coupling theory at high temperatures Q. Screening in other quantum number channels in 
the glue sector has also been studied [j]. Screening in colour singlet channels due to quark bilinear (meson- like) and 
trilinear (baryon-like) currents [l[ was understood as the first signal of deconfinement above the chiral symmetry 
restoring temperature in QCD with dynamical quarks |5]. Analyticity arguments relate these hadron- like screening 
masses in the low-temperature confined and chiral symmetry broken phase to the (pole) masses and properties of the 
hadrons. This has implications for models of heavy-ion collisions such as the hadron resonance gas model. 

One more application is to the viability of resummation of the weak coupling series at high temperature using 
dimensional reduction. This is possible only if the lowest screening mass belongs to the glue sector Q. It turns out 
that dimensional reduction does not work just above the chiral cross over temperature, T c . 

Hadron-like screening masses have been studied extensively 0,9- In QCD with light dynamical quarks they have 
been studied before using 2 flavours of staggered quarks @ and with 2+1 flavours of p4 improved quarks [loj]. They 
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have been studied also with overlap valence quarks and staggered sea quarks [TTJj . In all these studies the renormalized 
light quark masses are almost equal, and nearly physical. 

In this work, we extend previous studies through the analysis of meson-like spatial correlation functions in 2-fiavour 
QCD with staggered quarks. This brings the state of the art for dynamical staggered quarks into the regime of lattice 
spacings already reached using quenched overlap quarks. The organization of our paper is the following: In Section 
HTl we discuss operators selected for our analysis and the technical details of our fitting methodology. We investigate 
chiral symmetry restoration through the correlation functions in Section IIIII Our results on the screening spectrum 
are presented in Section [IVI along with a finite volume analysis. Finally we summarize the main details in Section fVl 
A technical point about the covariances of measurements of correlators is dealt with in Appendix [A] 



II. CONFIGURATIONS, MEASUREMENTS AND ANALYSIS 

A large part of this study uses decorrelatcd gauge configurations described in [13J. Two light flavours of staggered 
quarks were used with the bare quark mass tuned so as to give ~ 230 MeV at zero temperature. The lattice 
spacing was a = 1/(6T), i.e., N t = 6. The extraction of the temperature scale was explained in [l3j; we note that 
T c was identified there through the peak of the Polyakov loop susceptibility. The lattice volumes, V = (aN s ) 3 were 
set using 12 < N s < 24, i.e., the aspect ratio C = T\/V — N s /N t between 2 and 4. This provided basic control of 
finite volume effects; most of our results are reported for the largest volumes, £ = 4. In addition, wc performed a 
detailed finite volume scaling study at T — 0.94T C . For this we generated configurations with 4/3 < ( < 5, with all 
other parameters fixed as before. We measured autocorrelation times as before, and used at least 50 decorrelated 
configurations for our measurements. 

The meson screening correlation functions projected to zero momentum are — 

°1 = \ /l> [G(x,z)Gt(x,z)] 7 (x)^ (1) 

where x stands for sites labelled by the triplet (x, y, t), the number of terms in the sum is the same as the volume of 
such a slice, V = N x N y N t , G(x, z) is the inverse of the Dirac operator, i.e., the quark propagator from the origin to 
the point (x, z), the angular brackets denote an average over gauge field configurations with the correct weight, and 
the staggered phase factors 7 (x) pick out the quantum numbers, 7, of the meson under study. 

In this work we have taken all eight possible local staggered phases. At T = they would correspond to the flavour 
non-singlet scalar (S) (corresponding to the ao meson), the Goldstone pion (PS), and three components each of the 
local vector meson (V) and the axial vector ( AV) . Symmetry operations of the spatial slice interchange the components 
of the V and AV, so the three components are expected to be identical after averaging over gauge configurations. 

Since we measure spatial direction correlators at finite temperature, the symmetries of the (x, y, t) slice orthogonal 
to the direction of propagation are not the same as they would be in the corresponding zero temperature computation 
pj . The S/PS operators both lie in the trivial representation, called the A^, of the spatial direction transfer 
matrix. The sum of the x and y polarizations of the V/AV, and, separately, the t polarization, also lie in the A^ + 
representation. These six different kinds of Af + operators do not mix under the symmetries of the (x, y, t) slice, and 
hence we need separate notations for them. For the S/PS correlators it is economical to carry on the T — notation. 
For the sum of the x and y polarizations of the V we use the notation Vs (and AVs for the sum in the AV sector) 
and for the t polarizations we use the notation Vt and AVt. The difference of the x and y polarizations of the V/AV 
lie in a non-trivial representation called the Bf + . We use the notation VB and AVB for these. These particular 
realizations of the correlator have earlier been seen to vanis h fla| . 

We will also have occasion to use the S and PS susceptibilities [15[ defined as 

xps = J2 c * s > and xs = ^2(-iyc^ s . (2) 

z z 

The construction uses a fact that we demonstrate later: at high temperatures the S/PS correlators are essentially 
dominated by a single parity state. 

The inversion of the Dirac operator was done using a conjugate gradient (CG) algorithm, as usual. The tolerance 
for stopping was chosen such that the residual vector had squared norm less than eN t N 3 . We investigated whether 
we had an acceptable stopping criterion by monitoring 
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FIG. 1: S(e), see eq. @, as a function of the separation z for varying e, with a fixed value of e' = 10 



where C PS [e] is the PS correlation function obtained when the stopping tolerance parameter is e. We chose a fixed 
(' = 10~ 5 . In Figure Q] we show (5 computed on a randomly chosen test configuration at T = 0.94T C with £ = 4. The 
configuration to configuration variance of C is about 2-5% of the expectation value, so keeping S < 0.01 suffices. 
Clearly, the errors converge fast, and our choice of e = 10~ 5 is seen to be more than sufficient. 

Each staggered correlation function may contain contributions from two parity partners, and can be parametrized 
through the doubled-parity fit 

C(z) = Ai{e~^ z + e-" 1 ^*- z) ) + {-l) z A 2 {c-^ 2Z + e-M"*-*)^ (4) 

where /ii and y,i are the screening masses of the lightest natural parity meson appropriate for the operator used and 
its opposite parity partner. 

Since measurements of the correlation function at different distances, z, are made using the same gauge configu- 
rations, they are correlated, and the fit must take care of these correlations. Therefore, we used the definition of 

x 2 

zz' 

Here z is the spatial separation, C z are the measured expectation values of eq. ([l}, the function C(z) is the 2-mass 
form of eq. ((H), and Y, zz > is the covariance of C z and C z i, When T, zz r is diagonal, the definition reduces to the more 
familiar one. In actuality, the correlation coefficients are fairly high, so the matrix S 2Z / is nearly singular. The 
inversion was done in Mathematica to an accuracy of O(10~ 10 ). The errors in the inversion were therefore negligible 
compared to the statistical errors in the measurements, a z , which were of the order of a few percent. See Appendix 
[A] for further discussion of this procedure. 

A check on the consistency of the results obtained from fits is to use local masses. Due to the even-odd oscillations 
for staggered fermions, we used the definition of [9]- 

C z+1 = cosh[-m(z)(z + l- N z /2)} 

C z _i cosh[-m(»(z- 1- N z /2)]' [ ' 

Given the measurement on the left, the effective mass, m(z), can be extracted by solving the equation and errors 
estimated by jack-knife. This differs from a procedure where successive time slices are used for the modified correlator 

{-iyc z S3- 

In the chiral symmetry broken phase there is no particular relation between fix, [ii and A\, A 2 for different 
correlators. However, when chiral symmetry is restored, the staggered phases give 

d s = {-i) z c s z , c* Vs = (-i) z c^ s , c? vt = {-ifcY 1 . (7) 
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This implies the relations 

A\ s =A% Vs , fxY s =fi2 Vs and (Vs o AVs), (8) 

and similarly for the Vt and AVt or the S and PS channels. These relations are very easily demonstrated by using 
the projections 

C (±s) = c ps ± (_iy G f , Ci ±vs ^ = C v z s ± {~\) z Ct v \ Ci ±vt ^ = Cp ± {-l) z Ct vt . (9) 
If the correlators c| vanish for all z then chiral symmetry is restored for the full spectrum of excitations. 

III. THERMAL EFFECTS AND APPROXIMATE CHIRAL SYMMETRY RESTORATION 
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FIG. 2: The panel on the left shows Vs and Vt correlators at z — as a function of T/T c . Note the abrupt lifting of the 
degeneracy at T c . The panel on the right shows the ratio of the Vs and Vt correlators as a function of z, normalized by their 
values at z = 0. The data at T = 0.97T C is displaced slightly to the right for clarity. 



The screening correlators at any non-zero temperature should be decomposed according to the symmetry group 
of the finite temperature slice. At sufficiently low temperature, however, one expects the Vt and Vs correlators to 
be nearly equal, and the symmetries of the T — problem to be realized approximately. We investigated this by 
computing the ratios of the Vt and Vs correlators (normalized to be unity at z = 0). The statistical analysis was 
performed using a bootstrap, since the distribution of the ratio is not expected to be Gaussian The results below 
T c are shown in Figured] the ratio is consistent with unity at all z. The normalization is the value of the ratio of the 
correlators at z = 0. In Figure[2]we have plotted Cq s and C^ 4 as a function of T. Below T c the two are equal within 
statistical errors. The two facts taken together imply that the T = symmetries remain good until rather close to T c . 
Quite abruptly, just above T c this higher symmetry is broken, and the symmetry of the finite temperature problem is 
obtained. Similar results are obtained for the AVs and AVt. In view of this, in most of our subsequent analysis, we 
will group the correlators below T c into S, PS, V and AV. For T >T C we will continue to use the decomposition into 
S, PS, Vs, Vt, AVs and AVt. 

The next question which we examine is whether the correlation functions exhibit chiral symmetry restoration at 
any temperature. As discussed earlier, the most straightforward way to examine this is to plot d ^ at each T and 

ask whether it is consistent with zero at all z. In Figure [3] we show these quantities at two temperatures above T c . 

(-S) 

From the correlator Cs we see that at T — 1.48T C the symmetry is clearly restored, whereas for T = 1.012T C the 
symmetry is broken. The correlator c| Vs%> is consistent with zero for zT > 1, but at distances less than \/T there is 
clear chiral symmetry breaking close to T c . The correlator Vz Vt ^ most nearly exhibits chiral symmetry restoration 
immediately above T c , with only the value at z = being significantly non-zero. At higher temperatures all the C| 7 ^ 
are consistent with zero at all z, thereby indicating chiral symmetry restoration. 

In order to extend this analysis to all temperatures it is useful to introduce a less local quantity, 

zz' 
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FIG. 3: The chiral projections C, 



(-S) (j(-Vs) 



and Cz Vt ^ of eq. @ at two temperatures as computed on lattices with £ = 4. 



The Vt channel shows chiral symmetry is close to being restored immediately above T c . 
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TABLE I: This table lists the values of \ 2 a t different temperature for tests of the hypotheses that various correlators vanish. 
The number of degrees of freedom in all these cases is 12, since there are 13 independent values of z on the lattices with £ = 4 
with periodic boundary conditions. In order to rule out the hypothesis that a correlator vanishes at the 99% CL, the value of 
X 2 should be more than 36. 



where £ zz < is the covariance matrix of the measurements of the correlator at different distances. This is a measure of 
the likelihood that the correlators at all z are consistent with zero, and hence that chiral symmetry is restored. Values 
of these variables are collected in Table Q] Note that \vt shows a distinct change at T c , although it is consistent with 
chiral symmetry restoration only at T = 1.21T C . \Vs also shows a change at T c , although it is less dramatic. From 
Figure [3] it would appear that the change in Xvs a ^ is due to the long distance (z > 1/T) correlation function 
becoming consistent with zero, whereas the short distance (z < 1/T) part disappears only at larger T. Xs> 011 the 
other hand, does not seem to undergo any significant change at T c and signals chiral symmetry restoration only at 
T = 1.33T C . One sees the difference in behaviour in Figure [3] the S/PS correlators, unlike the V/AV, do not show 
any kind of effective long distance chiral symmetry restoration. This spatial structure has not been noticed before, 
and could be worth further investigation in future. 

The late restoration of chiral symmetry breaking can be understood from the fact that the non-vanishing quark 
mass provides explicit chiral symmetry breaking. In the chiral limit, there is a phase transition at T c . In the high 
temperature phase there is, effectively, a single scale, T, so the screening mass /i cx T. However, when there is a 
non- vanishing bare quark mass, m, there is no phase transition at T c but only a cross-over. In the absence of a phase 
transition, one could have fi/T — f{m 7I /T) where m T is the pion mass at T — (we have traded the bare quantity 
m for a renormalized measure of chiral symmetry breaking, m T ). At large T, when the argument of the function 
becomes small, / should go to a constant. On lowering T from large values, non-constant behaviour should become 
visible when the argument becomes of order unity, i.e., at T/T c ~ m^/T c . Since our simulations are performed with 
TOtt ~ 230 MeV and T c ~ 175 MeV, this argument implies that explicit chiral symmetry breaking should manifest 
itself up to T ~ 1.35T C , which is what we see. Such an argument would lead us to expect that in the real world chiral 
symmetry breaking in screening masses should not be visible above T c . This may have some bearing on the relation 
between T c defined through susceptibilities of the deconfinement and chiral order parameters (l7j . An alternative 
explanation of the late restoration of chiral symmetry is due to approximate restoration of the Ua(1) symmetry [To| . 
A future computation with different sea quark masses can easily distinguish between these two alternatives. 

The projection of the local V and AV channel correlators on the channel is expected to vanish. This was 

demonstrated in 0, [l6| with lattice spacing a — 1/ (4T) . Here we investigate the vanishing of these correlators at 
smaller lattice spacing using a correlated x 2 definition similar to that above. In Xvb/avb ^ ne f ac t° rs of Cj 7 in eq. 
(TTTJ1) are replaced by the VB or AVB correlator. The results are collected in Table U and show that the VB and AVB 
correlators vanish. 

In Figure @] we show the ratio of the PS and S susceptibilities (see eq. [2]). As expected, they become equal at 
T = 1.33T C , which is the point where the two correlators begin to satisfy eq. ©. At lower temperatures xps is larger, 
essentially because fips is smaller than fj,s- 

In the next sections we examine the other projection of the correlation function, c| +7 ' ) . This is non-zero at all 
temperatures. At high enough temperature one might expect the whole correlation function to be described in a 
weak-coupling theory. In Figure [5] we show the correlators at T = 1.92T C . One sees that the correlation function is far 

(-\-S) 

from the free field theory result, especially the correlator C z , indicating that even at this temperature the theory 
cannot be treated as weakly interacting. This is consistent with previous observations [9l-fli"j|. 
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FIG. 4: The ratio for susceptibilities of the pion to the scalar meson. 
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FIG. 5: The correlators Ci +7 ' obtained at T = 1.92T C with bare quark mass m/T = 0.1. The lines show the result for free 
field theory computed with the same quark mass. 

IV. SCREENING MASSES 

Figure |5] displays all four screening correlators at two temperatures, one each in the hadron and the plasma phase, 
i.e., at T = 0.97T C and T = 1.92T C . The V/AV correlators show clear even-odd oscillations at both temperatures, 
whereas these staggered artifacts are less clear in the S/PS correlators. This has clear implications for the fits: the 
former always requires a doubled-parity fit of the form given in eq. ((4]), whereas for the latter a single parity form 
may suffice. 

In Table Ql] we show that this is indeed correct. In the PS channel a single mass fit is acceptable judging by the 
value of x 2 /DOF, and the fitted value does not change significantly when a doubled-parity fit is performed. In fact, 
when a doubled-parity fit is attempted to the data, the mass of the parity partner is ill-determined. In the V channel, 
on the other hand, the doubled-parity fit turns out to be indispensable. In the AV channel the behaviour is similar to 
that in the V. In the S channel a single mass suffices, although one needs a (— l) 2 factor multiplying the exponential 
to take care of staggered oscillations. Only at the highest temperature were we able to extract a second mass from the 
S/PS correlators. The table also shows that the fitted masses are reasonably stable against changes in the fit range 
in both the PS and V channels. We find similar results for all T. 

Interestingly, although the covariance matrix is nearly singular (the smallest correlation coefficient being about 
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FIG. 6: The correlators for all the mesons at representative temperatures in the two phases on the lattice with £ = 4. Note the 
absence of even-odd oscillations for the S and PS correlators. 





PS 


V 






Uncorrelated single 


parity fit 




Range 


3-11 


4-11 


5-11 




4-9 


4-10 




A 


3.61(7) 


3.59(6) 


3.59(6) 










m 


0.401(2) 


0.400(2) 


0.400(2) 










X 2 /DOF 


0.24/7 


0.01/6 


0.01/5 












Correlated single parity fit 


A 


3.64(18) 


3.59(16) 


3.58(16) 




0.053(4) 


0.053(3) 




771 


0.401(3) 


0.400(3) 


0.400(2) 




0.484(7) 


0.480(8) 




X 2 /DOF 


11.0/7 


2.1/6 


1.9/5 




9.5/4 


9.8/5 






Correlated doubled-parity fit 
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A 1 




3.59(14) 


3.59(12) 


3.58(11) 


-0.67(4) 


-0.70(3) 


-0.55(5) 


Mi 
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1.20(3) 


1.14(3) 


A 2 




0.02(6) 


0.04(12) 


3.1(*) 


0.38(6) 


0.49(4) 


1.4(4) 


fJ-2 




0.9(*) 


0.9(*) 


1.5(*) 


1.62(8) 


1.71(5) 


1.98(11) 


x 2 /dof 




1.60/4 


1.58/3 


1.28/2 


6.55/6 


6.56/6 


2.20/3 



TABLE II: Fits to PS and V Correlators at T = 0.97T C on a lattice with £ = 4. An asterisk on a number indicates that the fit 
value is not determined reliably. 



0.8), the difference between the parameters extracted using or neglecting the covariance matrix in the PS channel is 
marginal. The major difference seems to be that the value of x 2 obtained when covariances are neglected are clearly 
too small for the usual statistical interpretation. 

In Figure [7] we demonstrate that the fitted masses agree with the local mass extracted from eq. ((6]) . We draw 
attention to the fact that the local masses exhibit a very well developed plateau, indicating that the correlation 
functions of fixed parity can be well described by a single mass. The spatial structure of chiral symmetry restoration 
shown in Figure [3] is visible also in the local masses at high temperature. 

In Table UTTl we collect the results for the fitted masses at all temperatures. We checked in all cases that the local 
masses for z > 1/T were compatible with these fits. For the V/AV correlators we also checked that if the fits were 
restricted to z < 1/T the fit results were generally different. Wherever doubled-parity fits are available, one can look 
for chiral symmetry restoration by checking whether or not the relations of eq. ([8]) are satisfied. Consistent with the 
analysis of Section llTll we find that this happens only for T > 1.33T C in the S/PS channels. Surprisingly, the equalities 
of eq. ([5]) hold in the V/AV channels, within statistical errors, from just above T c . 
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FIG. 7: The extracted local masses for the mesons at representative temperatures in the two phases on a lattice with f = 4. 
In the high temperature phase, where chiral symmetry is restored, only the PS, Vs and Vt masses are shown since the S, AVs 
and AVt exactly coincide. 
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FIG. 8: Screening masses as functions of T/T c . Below T c the S, PS, V and AV channels are shown. At T c and above the S, 
PS, Vt and AVt are shown. The free continuum value of 2tv is indicated with an arrow. 



In Figure [8] we plot the lowest screening mass in each channel as a function of T/T c . Above T c we could plot three 
channels. To avoid clutter we plotted only the S/PS and Vt/AVt channels. As one sees in Table Hill the lowest 
Vs/AVs masses are slightly larger, but consistent with Vt/AVt at the 2-a level. All the features discussed are clearly 
visible here. Also visible is the fact that [ips/T increases monotonically with T whereas Hs/T dips near T c . Note 
also that /^y t /T may approach its ideal gas value from above, becoming consistent with the limit already at T ~ 2T C . 
However ns/T remains about 20% below this limit even at the highest temperature we explored. We shall return to 
this point later when we discuss the continuum limit. 
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TABLE III: Screening masses at different temperatures. The fit range was z/a=3-ll except for the S/PS at 1.92T C , where 
the fit range was z/a=4-ll (the larger range gave very large \ 2 without appreciably changing the best fit values). Also note 
that except at 1.92T C , the PS fit was done with a single-parity fit form, and therefore has two more degrees of freedom than 
the other channels. A dash indicates that some mass could not be obtained because staggered oscillations were not visible. At 
temperatures below T c the analysis was performed on the V/AV channels. 



A. The role of explicit chiral symmetry breaking 

In free field theory one has no pion and the explicit chiral symmetry breaking scale is the quark mass. In this case 
one has 



T 



2tt 



1 



1 / m \ ; 

2 \rfT) 



(11) 



Substituting the bare quark mass into this expression, it can be seen that the effect is of the order of a few parts in 
10 5 , and hence negligible. However, it turns out that in weak-coupling theory one has to insert in the above equation 
the thermal mass of the quark Since this is gT/y/3, and g is large near T c [lj|, the effect can be significant. Of 
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FIG. 9: The valence quark mass dependence of various quantities at T = 1.21T C when the bare sea quark mass is held fixed at 
the value m/T c = 0.1 (corresponding to ~ 230 MeV). The screening masses in the S/PS sector depend significantly on the 
valence quark mass, whereas the V sector screening correlators are insensitive. 
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FIG. 10: Local masses in the PS channel for different valence quark masses at T = 1.33T C with £ = 4. The horizontal lines are 
the 1-sigma bounds from the fit to the correlator for m va i = m sea . 



course, when g is large, the weak-coupling theory is unlikely to be quantitatively useful, and should be taken only as 
an indication that one must explore the quark mass dependence of the chiral symmetry breaking seen in the screening 
masses. 

Changing the quark mass involves in principle a completely new set of computations, and lies beyond the scope of 
this paper. However, it is possible to change the valence quark mass (m va i) without changing the sea quark mass 
(m S ea)- As a result, such a restricted study can be done without a large cost in CPU time. In view of this, we repeated 
our analysis above T c with m va i chosen to be 2/3 and 1/3 of m sea . Sample results are shown in Figure |H1 

We find that a change in the valence quark mass has insignificant effect in the V/AV channels (see, for example, 
the second panel of Figure [9]) . However, there are statistically significant changes in the S/PS sector. Both the S and 
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PS screening mass increase with m va i. The difference aiso increases, although the limiting value for m va i = is finite. 
We find that 

= 0.46 ± 0.03 (for m sea /T c = 0.1, m val = 0). (12) 

T/T c =1.21 

It is not inconceivable that decreasing m sea in a future computation can lower this mass difference further. Note, 
however, that /ips/T and ns/T both drop as the quark mass changes, and move further away from the weak-coupling 
expectation. 

In fact, the dependence of correlation functions on the bare quark mass can be deeper than what was discussed 
above. We illustrate this in Figure [TUl where local masses in the PS channel are shown for several values of m va i- It 
seems that with decreasing m va i the plateau in the local masses becomes less well developed; for the smallest m va i, 
in fact, a distinct slope is visible. 

B. The effect of finite volume 
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TABLE IV: Screening mass estimates at T = 0.94T C at fixed a — 1/(6T) with changing spatial volume (£/T) 3 . 

Finite volume effects in the high temperature phase of QCD have been explored earlier Q and are by now well 
understood. These effects are under control as long as jips is large compared to the inverse size of the smallest lattice 
dimension, i.e., T. Since we find [ips/T > 2.5 in the high temperature phase, we expect that finite volume effects are 
under control for T > T c . 

In this paper we report on the magnitude of finite volume effects, below T c . We studied the screening correlators 
on lattices with 4/3 < £ < 5 at T = 0.94T C . As shown in Table ITVl and Figure fTTT finite volume effects are invisible 
within statistical errors. Again, since l^ps/T ~ 2.3, this may not seem unexpected. 

It is worth the remark that while such studies can have some bearing on decay widths at finite temperature, much 
larger lattices and smaller quark masses and lattice spacings may be required for those. For example, the scalar 
under study cannot decay into two pions, and must have at least three pions in the final state. This is ruled out 
kincmatically at present. 

On the other hand, since ms/mps — 2-2, one must ask whether the long-distance behaviour of the S correlator is 
determined by a single scalar exchange or mult i- particle exchanges (this is the finite temperature analogue of particle 
decays, and we shall save space by using the word decay). In the continuum theory, this non-isosinglet scalar cannot 
decay into two pions. From this point of view, since ms/mps — 2.2 < 3, one could expect that the scalar does 
not decay. However, with staggered quarks on a lattice, there are spurious two pion states (taste multiplets) through 
which the scalar current could be correlated [2(| • The featureless behaviour of the volume dependence of the screening 
masses implies that no decays occur. However, this is a weak statement, because a 7% taste symmetry breaking is 
sufficient to forbid this decay on the largest lattice which we used. 

C. The continuum limit 

In Table fVl we extract the values of the screening masses in units of temperature, [ijT for two different temperatures 
in the chirally symmetric phase, at two different lattice spacings. From two pieces of data in each case it is hard 



MS ~ MPS 
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FIG. 11: Screening masses in various channels as a function of £ at T = 0.94T C . Note that there is nothing special about the S 
channel masses. 
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TABLE V: The lattice spacing dependence of fi/T in various quantum number channels. The temperature scale is rounded off. 
Data for a = 1/(4T) is taken from 0. 

to construct a continuum extrapolation. However, we can test whether the extrapolation is consistent with the 
expectation /i ~ 2irT by attempting a fit to the form 



= 2tt- 



N t 



(13) 



It is possible to get reasonable fits in the V/AV channel, yielding 



-13 ±2 (T=1.5T C 
-9.0 ±0.2 (T = 2T C ) 



(14) 



This is consistent with previous results. However, in the S/PS case this procedure fails to yield the expected result; 
the extrapolated screening length remains below the ideal gas value. One cannot rule out the possibility that the 
weak-coupling result emerges at even smaller lattice spacings. 

We note however, that the screening of meson-like correlations in the weak coupling theory occurs, not through the 
exchange of a single particle, but through multiparticle exchange. As a result, the zero-momentum correlator is not 
expected to be strictly exponential, but to have some curvature. Such a curvature was actually seen in computations 
in the quenched theory using staggered [2l[ and Wilson Q quarks using much smaller lattice spacings. An easily 
observed feature of such a curvature is that local masses do not show a plateau, but change continuously with z. Such 
behaviour was neither seen here (see Figure [7]) , nor with p4 improved quarks at the same renormalized quark mass 
(loj . However, as pointed out in Section TlV A[ when the quark mass is lowered such a feature could emerge. 
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V. SUMMARY 

We studied screening correlation functions and screening masses of meson-like probes of strongly interacting matter, 
both in the low temperature hadron and high temperature plasma phases. We used configurations described in [l3j ; 
these were generated using two flavours of dynamical staggered quarks with masses tuned to give m T ~ 230 MeV. 
Most of the results come from lattice spacing a = 1/ (6T) , although we have attempted to check assumptions about the 
continuum limit using earlier measurements with the same renormalized quark mass and lattice spacing a — 1/(4T). 
We have explicitly checked for finite volume effects, and found that these are negligible when we take the aspect ratio 
C = 4. 

We checked that the correlators at T < T c effectively have the symmetries of the T = transfer matrix, and that 
there is a fairly abrupt transition at T c to the symmetries of the screening transfer matrix (see Figure [2]). Although 
the QCD cross over occurs at T c , we found a lack of parity doubling in the spectrum of screening masses up to a 
temperature of 1.33T C (see, for example, Figure 0] and Table [I]). We argued that explicit chiral symmetry breaking 
due to the quark mass is visible up to a temperature of 1.35T C , as seen here. This argument leads us to believe that 
in the limit of realistic quark masses, such a breaking would not be visible above T c . Interestingly, there turns out to 
be a lot of structure in this apparent breaking of chiral symmetry above T c . The V/AV correlators are equal to each 
other at distances z > 1/T, and the chiral symmetry breaking in this sector is entirely a short distance effect (see 
Figure [3]) . In the S /PS sector the non-degeneracy of the correlators persists into the long distance regime. 

Non-degeneracy of the S/PS correlators could also be due to £/a(1) symmetry breaking. This is suspected to persist 
well into the high temperature plasma phase [HI . We tested what happens to the S/PS difference as the valence 
quark mass is changed (see Figure [5]). Our results imply that other physics effects can be disentangled from the 
explicit symmetry breaking effect due to finite quark mass only through computations with smaller quark masses. 
Alternatively, one could examine the short distance part of the V/AV correlators for signals of such microscopic 
physics. 

We made the first study of hadron decays at finite temperature (below T c ) through a systematic exploration of the 
volume dependence of screening masses. We found no significant volume dependence (see, for example, Table ITVl and 
Figure [TT]), indicating the stability of the scalar. As we discussed already, this study needs to be carried out with 
smaller quark masses so that us/ Hps > 3, or at smaller lattice spacings, so that taste violations are reduced. 

We combined the analysis of this paper with data from an earlier source Q in which the same renormalized quark 
masses were used to study screening at a coarser lattice spacing, a = 1/(4T) to explore the continuum limit. Being 
restricted to only two values of the lattice spacing at each T, we ask whether the continuum limit of screening masses 
is compatible with the ideal gas expectation, 2ttT, in the high-temperature phase. We find that it is, in the V/AV 
channels, but not in the S/PS channels. 

If the high temperature phase is deconfined, then correlations of static currents with meson quantum numbers must 
be mediated by the exchange of a quark anti-quark pair. The most straightforward signal of this is that the local 
masses do not show a well-developed plateau. In most of our studies we did not see this. Only in a study with rather 
small valence quark masses did we see a signal of such behaviour (see Figure I10[) . Studies with lower sea quark masses 
in the future will be needed to resolve the question of deconfinement above T c in QCD with physical quark masses. 
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Appendix A: Covariance matrices and statistical data compression 

The correlation functions are evaluated on a fixed set of gauge configurations. In general, the estimates of the 
correlator at different separations are not statistically independent. The covariance matrix of the correlators measures 
the degree of independence, and is defined as 

£ 2Z , = <(C? - <C?» (C], - {CI))) , (Al) 

where the averages are over the ensemble of gauge configurations. The diagonal elements are the variances. The 
matrix is positive and symmetric by construction, so that it has well defined (positive) eigenvalues and orthogonal 
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FIG. 12: The measured correlator (red circles) and the compressed version obtained by projection on to a single eigenvector of 
E (blue squares) for the PS and V at T = 1.92T C and 0.97T C . 



eigenvectors. It is found that it often has large condition number, i.e., the ratio of the largest and smallest eigenvalues 
is large. Any fits to the correlation function are made by minimizing the x 2 -function 

X 2 = D*YT X D, where D = C - f, (A2) 

C 1 is the vector of measured correlation functions, with components C] , and / is the corresponding vector made 
from the function to be fitted. 

The covariance matrix is often inverted using a singular value decomposition, E = UDV T where U and V are 
orthogonal matrices (in this case U — V) and D is diagonal. The inverse matrix is usually obtained as E -1 = VD^U 7 ', 
with the largest components of D _1 set to zero [231 ] - This procedure is used to prevent errors arising from finite 
precision rounding. For this part of the analysis we use Mathematica, which allows arbitrary precision computations. 
By tuning the precision of the computation from 10 digits to 40 digits, we checked for control of rounding errors. Our 
conclusion is that the rounding errors can be kept under control without setting large numbers to zero. As a result, 
in our ensuing investigations we are assured that the eigenvalues of E are the inverses of those of E -1 . 

The covariance matrix and its eigenvectors are determined essentially by the statistical properties of the measure- 
ment and not directly by the correlation function. One extreme can be illustrated by the imagined case where we 
generate a very large set of decorrelated gauge configurations, and use a disjoint subset for evaluation of the correlator 
at each z. Since the measurements at each of the z are then statistically independent of each other, E must be 
diagonal. As a result, the eigenvectors have support only on a single separation z. 

Take the eigenvalues of E _1 to be X a with corresponding eigenvectors v a . Each correlation function can be written 
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as a linear combination of the v a , i.e., 

C 7 (A3) 

a 

Our purpose in doing this is that if some of the c 7 turn out to be small then one can perform a noise reduction by 
dropping the small terms. Since c 7 = C 7 • v a , the errors on the dot product are clearly given by 

cr 2 (c 7 ) = v az v azl Yj zz , = — . (A4) 

In writing this formula we have taken v az to have no errors. This approximation can be removed, if necessary, by a 
bootstrap analysis of the errors in c 7 . If many of the components, c 7 , are zero within errors, then we can drop them 
from the analysis. In this case we say that the correlation function is compressible. Such a notion corresponds to that 
of 'lossy data compression' in a variety of contexts, including jpeg image compression and mp3 audio compression. 

In the imagined case which we discussed before, generically all the c 7 will be non-zero and the correlator will not 
be compressible. Only if some of the measurements are so noisy that they are compatible with zero can we drop them 
from the fits. This example illustrates the connection of compressibility with that of statistical independence. 

The actual situation we investigated was the usual one where the correlation function at all separations were 
measured on all configurations and the lowest mass was obtained by a fit to the long-distance correlation function 
using the cuts mentioned in Section ITV1 i.e., 3 < z/a < 13. Very surprisingly, we found that the PS correlator can 
be compressed to exactly one component, corresponding to the smallest X a both above and below T c . This means 
that the corresponding eigenvector v a contains the full information on the correlator. The V correlator can also be 
compressed down to a very small number of components, usually one or two. These surprising results are shown in 
Figure [121 When we examine the correlation functions at all distances, then again we find a high, but lesser, degree 
of compressibility with two or three eigenvectors being needed for the description of the data. 

What is the origin of these extremely strong covariances in the data? The gauge configurations we used were 
chosen so that thermodynamic operators such as the action, Polyakov loops, chiral condensate, various quark number 
susceptibilities, are decorrelated between one configuration and another. Several of these quantities, for example, 
those which involve fermions, are highly non-local. In view of this, autocorrelations between configurations are not 
the cause of the high compressibility of the data. 

So it must be the correlation function itself which generates these covariances, with the result that one of the 
eigenvectors, v a contains all the information present in the measured correlator. Could this be due to the fact that 
the Dirac propagators are highly non-local? If so, the baryon correlators must also be highly compressible, whereas 
glueball correlators need not be. Or could it be that the relatively low masses are at the root of the covariance? If 
so, neither the baryon nor glueball correlators need be compressible. Since correlated fits are standard technology 
in the fitting of masses, measuring statistical compressibility of correlators is a trivial extension of standard lattice 
measurements. It would be interesting to have more data in future on the compression of correlation functions. 
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